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ABSTRACT 


Structures  are  calculated  for  Magnesium  Phthalocyanine  (MgPc)  and  its  radical  anion 
doublet  (MgPc-),  using  both  ab-initio  (6-3 1G**)  and  semiemperical  (INDO/1)  self-consistent 
field  approaches.  The  anion  is  first  order  Jahn-Teller  distorted,  and  the  various  distortions  that 
are  possible  are  examined.  The  electronic  absorption  spectra  of  both  molecular  species,  and  the 
effect  varying  the  degree  of  distortion  has  on  the  computed  anion  spectrum  are  discussed.  These 
results  suggest  that  the  four-orbital  model  often  applied  to  porphyrin  systems  in  interpreting  the 
low  energy  spectrum  is  incomplete  for  the  anion  case.  We  further  conclude  that  the  Jahn-Teller 
distortions  calculated  by  either  ab-initio  or  semiemperical  models  may  be  too  great. 
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Introduction 

Recent  articles1"*  concerning  the  experimental  UV-VIS  spectral  properties  of  the  radical 
anion  of  metal  phthalocyanines  and  their  interpretation  have  spurred  our  interest  in  the  calculation 
of  such  properties.  In  this  work  we  obtain  calculated  UV-VIS  absorption  spectra  for  both 
magnesium  phthalocyanine  (MgPc)  and  its  anion  radical.  Workers  in  this  area  have  explained 
the  spectral  features  of  the  phthalocyanines  using  the  four-orbital  model  of  Gouterman.5  It 
seems  reasonable  that  this  model  can  be  applied  to  the  electrically  neutral  phthalocyanines  as  the 
principal  structural  features  present  in  porphyrin,  the  species  for  which  this  model  was  developed, 
are  also  present  in  phthalocyanine.  Phthalocyanine,  like  porphyrin,  can  be  likened  to  sixteen 
member  conjugated  macrocycles  and,  usually,  possesses  a  two-  or  four-fold  principal  symmetry 
axis.  The  four-orbital  model  does,  indeed,  explain  many  of  the  features  of  the  low  energy  ix  ->  7r* 
transitions  localized  to  the  porphyrin-like  macrocycle.  It  is  not  obvious  that  this  elegantly  simple 
model  can  be  applied  to  anionic  structures,  where  the  reducing  electron  occupies  the  lowest 
energy  7r*  orbital  of  the  neutral  species.  We  will  address  the  applicability  of  this  model  to  these 
systems  through  the  interpretation  of  our  calculated  spectra.  We  also  compare  our  calculated 
results  to  those  obtained  experimentally  in  an  attempt  to  better  understand  the  observed  spectrum. 

Before  a  spectrum  can  be  calculated  a  “reasonably”  good  structure  must  be  obtained  either 
by  direct  experimental  observation  or,  when  such  information  is  not  available,  by  computational 
techniques.  In  this  paper  we  will  compare  the  computed  spectra  of  four  quantum  chemical 
optimized  structures  to  spectra  acquired  experimentally,  as  the  calculated  spectrum,  especially  of 
the  anion,  is  remarkably  sensitive  to  the  geometry.  From  our  results  we  comment  on  the  ability 
of  theory  to  produce  good  structures,  and  to  predict  the  UV-VIS  spectral  properties. 

Computational  Methods 

Structures.  The  initial  MgPc  structure  was  obtained  through  the  use  of  a  Tektronics  Inc. 
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CAChe  workstation.6  The  neutral  structure  was  generated  and  then  refined  using  a  modified 
MM2  force  field7  available  with  the  CAChe  system.  This  classically  optimized  structure  then 
served  as  the  starting  point  for  the  quantum  chemical  structural  optimization. 

For  the  ab-initio  structural  optimizations  we  used  the  GAMESS8  program  package.  We  opti¬ 
mized  both  the  neutral  and  the  radical  anion  at  the  Restricted  Hartree-Fock9  (RHF)  and  Restricted 
Open-Shell  Hartree-Fock10-15  (ROHF)  Self-Consistent  Field16  (SCF)  level,  respectively,  using 
a  6-3 1G**  basis  set.17-20  At  the  semiempirical  level  of  approximation  we  used  the  ZINDO21-30 
program  package  to  obtain  the  optimized  structures.  Once  again  both  structures  were  optimized 
at  the  SCF  level  using  RHF  for  the  neutral  and  ROHF31  for  the  doublet  anion.  These  SCF 
optimizations  were  carried  out  at  the  Intermediate  Neglect  of  Differential  Overlap/ 1 32  (INDO/1) 
level  of  approximation. 

In  both  cases,  ab-initio  and  semiempirical  models,  it  was  assumed  that  the  neutral  structure 
had  D4h  symmetry.  Adding  an  electron  to  the  neutral  D4h  species,  to  form  the  initial  structure  of 
the  radical  anion,  leads  to  the  single  unpaired  electron  populating  a  degenerate  eg( 7r*)  molecular 
orbital.  We  then  optimized  the  anion  radical  in  the  following  two  ways;  (1)  one  electron  in  one  of 
the  two  eg (tt*)  orbitals,  allowing  the  molecule  to  Jahn-Teller33  distort  into  two  possible  structures 
both  with  D2h  symmetry,  and  (2)  one  electron  averaged  over  the  two  eg(- k*)  orbitals.  This  latter 
averaging  method  maintains  the  D4h  symmetry  of  the  neutral  species  while  relaxing  the  structure. 
We  take  the  difference  in  energies  between  (1)  and  (2)  to  be  the  Jahn-Teller  distortion  energy. 

Electronic  Spectra.  Having  obtained  the  optimized  MgPc  and  MgPc-  structures  via  the 
above  described  methods,  we  then  utilized  them  to  obtain  the  calculated  UV-VIS  absorption  spec¬ 
tra,  using  the  ZINDO  program  package.  The  spectroscopic  INDO  Hamiltonian  (INDO/S)21, 22 
was  used  to  generate  a  set  of  molecular  orbitals  (MO)  which  were  then  used  to  construct  the 
multi-electron  basis  for  a  Singles-Only  Configuration  Interaction  (CIS)  calculation34  to  deter¬ 
mine  the  state  energies  and  absorption  intensities.  The  INDO/S  method  has  been  specifically 
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parameterized  to  yield  good  spectroscopic  transition  energies  at  the  CIS  level. 

For  the  neutral  species  the  active  orbitals  in  the  Cl  were  the  50  lowest  energy  virtual 
molecular  orbitals  and  the  50  highest  energy  occupied  molecular  orbitals.  This  same  active 
space  was  spanned  for  the  CIS  calculation  of  the  anion.  The  active  space  may  be  adjusted  up  or 
down  by  one  orbital  as  needed  to  prevent  the  splitting  of  a  degenerate  MO  level,  or  in  the  case  of 
the  D2h  structures,  those  orbitals  that  would  be  degenerate  if  the  structure  were  of  D4h  symmetry. 


Results 


The  atomic  numbering  scheme  used  to  label  the  symmetry  unique  bond  lengths  and  bond 
angles  is  given  in  Figure  1 .  Their  calculated  values,  for  both  optimization  methods  and  both 
molecular  species,  can  be  found  in  Tables  I  and  II. 

[PLACE  FIGURE  1,  TABLE  I,  AND  TABLE  E  HERE.] 

The  calculated  UV-VIS  absorption  spectra  for  the  neutral  species  are  tabulated  in  Table  V 
for  the  INDO/1  optimized  structure  and  in  Table  VI  for  the  ab-initio  SCF  optimized  structure; 
the  calculated  stick  spectra  were  also  fitted  to  Lorentzian  distribution  functions  and  plotted.  (See 
the  Appendix  for  a  complete  description  of  the  Lorentzian  function  used.)  These  plots  can  be 
found  as  Figures  2 a  and  2b.  For  the  radical  anion  the  calculated  UV-VIS  absorption  spectrum 
of  the  INDO/1  structure  can  be  found  in  Table  VE  and  the  results  for  the  ab-initio  SCF  structure 
can  be  found  in  Table  VIE.  Their  Lorentzian  distributions  are  plotted  in  Figures  5a  and  5b. 

Those  molecular  orbitals  which  participate  in  the  description  of  the  electronic  excitation 
spectra  are  ordered  by  energy  and  labeled  by  symmetry  in  Figure  3a  for  the  INDO/1  optimized 
neutral  structure,  and  in  Figure  3 b  the  anionic  ENDO/1  structure  is  likewise  labeled.  Finally, 
Tables  I  and  E,  and  Figure  1  list  and  show  how  the  optimized  anionic  structures  Jahn-Teller 
distort. 
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Discussion 


Neutral  MgPc  (D4h).  The  vapor-phase  absorption  spectrum  of  neutral  MgPc  has  been 
reported  many  times.35, 36  The  five  distinct  peaks  between  800  nm  and  200  nm  have  been 
labeled  as  Q,  B,  N,  L  and  C,  the  same  description  used  for  porphyrin  spectra.35  We  adhere 
to  this  convention  while  explaining  our  results.  The  experimentally  observed  bands,  as  well 
as  those  calculated,  can  be  found  in  Table  IE.  Our  assignment  of  Q,  Bi,  B2,  N,  L  and  C  is 
somewhat  arbitrary.  We  assign  these  labels  to  the  six  lowest  energy  Eu  bands  of  observable 
intensity  for  both  the  INDO/1  structure,  and  the  ab-initio  structure.  For  both  calculated  structures 
the  band  labeled  B  is  composed  of  two  calculated  Eu  transitions.  For  both  structures  these  two 
bands  lie  too  close  in  energy  to  be  readily  distinguishable  experimentally,  see  Figures  2  and 
6a;  in  the  literature  these  have  been  noted,  and  labeled,  as  Bi  and  B2.  The  split  B  band  is 
not  resolved  in  the  gas-phase  spectrum  of  Reference  35,  but  is  seen  in  the  MCD  spectrum 
of  Reference  4.  For  the  most  part  the  Eu  bands  are  described  by  excitations  localized  to  the 
macrocycle.  It  is  then  interesting  to  note  the  differences  in  the  calculated  Eu  bands  between  the 
two  optimized  structures,  and  correlate  these  observations  with  the  main  structural  difference. 
Although  there  are  several,  the  principal  ones  are  the  Mg  to  N  bond  distance,  2.098  A  vs  1.993 
A,  and  the  circumference  around  the  macrocycle,  21.872  A  vs  21.344  A  for  the  INDO/1  and  ab- 
initio  structures  respectively.  Although  we  have  not  been  able  to  locate  an  experimental  Mg-N 
bond  length  in  MgPc,  Mg-N  bond  lengths  in  porphyrins  and  chlorins  are  generally  found  to  be 
between  2.05  A  to  2.13  A.  The  effect  of  this  upon  the  calculated  spectrum  is  observed  in  Figure 
2  and  in  Tables  V  and  VI.  We  observe  that  the  calculated  excitation  energies  remain  relatively 
constant  between  the  two  structures  used  in  their  computation,  although  there  is  a  small  red 
shift  for  the  larger  INDO/1  structure,  while  the  computed  oscillator  strengths  differ  significantly; 
that  is,  the  band  centers  appear  insensitive  to  small  structural  changes,  while  oscillator  strengths 
tend  to  be  more  sensitive. 
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As  can  be  seen  in  Table  HI  there  is  some  discrepancy  between  the  reported  experimental 
spectra;  this  difference  is  larger  than  might  easily  be  explained  from  solvent  shifts.  The  positions 
of  the  bands  of  the  solvated  spectrum  where  obtained  via  deconvolution  calculations.35  Our 
quantum  chemical  calculations  are  best  compared  with  the  observed  vapor-phase  spectrum. 

[  PLACE  TABLES  HI,  IV,  V  AND  VI  HERE  ] 

Most  of  the  intensity  in  both  porphyrins  and  phthalocyanines  derives  from  the  four-orbital 
model,  leading  to  four  excitations:  laiu  — *■  1  e*(xz),  1  a2u  — ►  1  e*(xz),  la\u  — >  1  e*(yz) 
and  la2u  —*■  le*(yz),  see  Figure  3.  For  mg-porphyrin  the  lai„  and  la2U  molecular  orbitals 
are  calculated  to  be  nearly  degenerate  as  are  the  excitations  from  these  orbitals  into  the  le*. 
These  strongly  allowed  excitations  mix  through  the  configuration  interaction  treatment.  For  the 
lower  energy  transition,  Q  band,  the  intensities  cancel,  and  for  the  higher  energy  B  band,  they 
reinforce.  In  the  magnesium  porphyrin  INDO/S-CIS  calculation  the  Q  band37  is  composed  of 
60%  la\u  — ►  le*  and  40%  la.2u  —■ ►  le*,  and  the  B  band  is  composed  of  the  same  two  one- 
electron  orbital  excitations  with  the  ratios  reversed.  (We  number  the  molecular  orbitals  as  in 
Reference  36.  i.e.  The  HOMO  and  LUMO  are  given  the  label  1,  then  the  occupied  orbitals 
are  incremented  as  they  decrease  in  energy  and  the  virtual  orbitals  are  incremented  as  they 
ascend  in  energy.)  Calculations  with  higher  excitations  (doubles  and  triples)  contribute  to  the  B 
band,  reducing  its  calculated  transition  energy  and  oscillator  strength,  but  the  general  four-orbital 
picture  is  still  generally  valid,  especially  in  its  ability  to  describe  the  intensities.  Table  IV  shows 
that  for  both  Pc  structures  the  band  labeled  Q  is  comprised  of  approximately  90%  laiu  — ►  le* 
and  6%  la2u  -+  le*.  This  unequal  mixing  can  be  explained  now  by  the  greater  lai„  -  la2u 
orbital  splitting  found  in  the  Pc  calculation  of  0.097  eV,  compared  to  the  splitting  found  in 
the  MgP  complex  of  0.016  eV.  This  leads  to  strongly  allowed  excitations  into  the  le*  orbitals 
that  are  not  as  degenerate  as  in  the  porphyrin  case,  and  a  greatly  reduced  configuration  mixing. 
The  transition  intensities  do  not  cancel  for  the  lower  energy  band,  leading  to  Q  and  B  bands 
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of  near  equal  oscillator  strength,  as  is  observed.36  The  phthalocyanine  B2  band  is  composed 
of  63%  la2u  ->  le*  for  the  INDO/1  structure,  while  83%  of  this  orbital  excitation  is  shared 
between  the  Bj,  B2  and  N  bands  for  the  ab-initio  structure.  The  four-orbital  model  still  explains 
the  origins  of  the  Pc  Q  and  B  bands,  but  the  strong  mixing  between  the  1  a\u  ->■  le*  and 
la2u  le*  orbital  excitations  found  in  metal  porphyrin  is  not  present  in  our  MgPc  calculation. 

As  can  be  seen  from  Tables  V  and  VI  the  Q  band,  when  calculated  with  either  structure, 
is  approximately  1300  cm'1  too  low  in  energy  when  compared  to  experiment.  The  B  band  is  a 
different  story;  both  Ough36  et  al  and  Stillman4  et  al  report  two  close  lying  bands  in  the  region 
of  their  spectra  labeled  as  the  B  band.  Our  ab-initio  structure  places  two  strong  Eu  bands  within 
300  cm' 1  (4  nm )  of  each  other,  see  Table  VI,  centered  at  32900  cm'1  (304  nm).  A  third  strong 
Eu  band  lies  1500  cm"1  above  this  band  pair  and  when  plotted  in  Figure  2 a  forms  a  broad  peak. 
We  assign  this  peak  to  the  observed  B  band.  The  INDO/1  structure  places  two  strong  Eu  bands 
within  1300  cm f1  (14  nm)  of  each  other  centered  at  31500  cm f1  (318  nm).  When  plotted  this 
band  pair  also  forms  a  broad  peak,  see  Figure  2b,  and  we  assign  it  to  the  observed  B  band. 

It  is  of  interest  to  note  how  the  calculated  configuration  mixing,  and  thus  intensity,  vary  as 
a  function  of  small  structural  differences.  The  summed  oscillator  strengths  for  the  Q-N  bands  is 
6.4  for  the  ab-initio  structure  and  6.7  for  the  INDO/1  structure,  or,  approximately  the  same  for 
both  structures.  Inspection  of  Tables  V  and  VI  and  Figure  2  shows  how  the  intensity  is  distributed 
between  the  states.  Ten  Eu  bands  lie  below  45200  cm'1  for  both  structures.  Guided  by  their 
proximity  and  intensity,  we  have  assigned  the  calculated  bands  to  those  observed.  The  absorption 
envelopes  of  Figures  2 a  and  2b  both  contain  all  the  salient  features  of  Gouterman’s35  vapor-phase 
MgPc  spectrum,  reproduced  as  Figure  6a.38  In  our  opinion,  based  on  the  experimental  results 
of  References  35  and  36,  the  calculated  spectrum  of  the  INDO/1  structure  best  reproduces  the 
observed  MgPc  spectrum.  We  feel  this  is  supported  by  the  number  of  bands  with  appreciable 
intensity  and  their  proximity  to  those  observed  experimentally,  see  Table  HI. 
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It  is  of  interest  to  note  that  one  state  is  calculated  to  have  appreciable  charge  transfer  character 
below  50000  cm'1,  it  is  forbidden  and  of  A2g  symmetry.  For  both  structures  this  state  falls  at 
approximately  32000  cm'1  and  is  of  50%  Pc(tt)  -*•  Mg(pz)  character.  If  populated,  this  state 
should  lead  to  a  large  Mg  out  of  plane  displacement  as  Mg+  is  a  much  larger  ion  than  Mg+2.39, 40 

Jahn-Teller  distortions.  The  LUMO  of  neutral  MgPc  is  of  eg(7i*)  symmetry.  Placing 
the  reducing  electron  into  the  LUMO  of  the  D4h  structure  will  cause  the  structure  to  undergo 
a  first-order  Jahn-Teller  (JT)  distortion33  along  one  of  the  two  degenerate  vibrational  modes 
transforming  as  eg.  (We  define  the  four-fold  axis  to  be  the  z-axis,  then  the  plane  of  the  molecule 
defines  the  xy  plane.)  In  the  D4h  point  group  eg  vibrations  transform  as  do  xz  and  yz.  Thus 
one  expects  the  structure  to  distort  along  the  x-  and  y-axes  assuming  no  net  out  of  plane  forces, 
which  we  checked  for  and  did  not  find.  There  are  two  unique  choices  for  the  x-  and  y-axes,  one 
choice  has  the  aza  nitrogens  defining  these  axes,  and  the  second  choice  is  defined  by  the  pyrrole 
nitrogens.  This  situation  is  schematically  shown  in  Figure  4. 

If  we  choose  the  pyrrole  nitorgens  to  define  the  x-  and  y-axes  and  proceed  to  optimize  the 
structure  via  the  INDO/1  or  ab-initio  methods  we  obtain  a  structure  which  is  a  true  minimum 
on  the  electronic  potential  energy  surface.  While  we  did  not  have  the  computational  resources 
to  compute  the  Hessian  for  such  a  large  structure,  we  did  though,  in  various  ways,  distort  it 
from  the  computed  “equilibrium”  structure  and  in  all  cases  upon  reoptimization  the  structure 
returned  to  the  equilibrium  one.  In  this  case  the  D4h  structure  distorted  with  two  opposing  fused 
pyrrole-benzene  groups  moving  closer  to  the  central  Magnesium,  and  the  perpendicular  pair  of 
fused  pyrrole-benzene  groups  moving  away,  with  the  aza  nitrogens  now  forming  the  vertices  of 
a  square,  i.e.  the  diagonals  formed  by  connecting  N3  to  N5  and  N2  to  N4  are  of  equal  length 
and  are  perpendicular  to  one  another,  see  Table  I  and  Figure  1.  It  was  with  these  D2h  structures 
that  we  proceeded  to  calculate  the  UV-VIS  spectrum  of  the  anion. 

If  we  let  the  aza  nitrogens  define  the  x-  and  y-axes  and  proceed  to  optimize  the  structure 
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via  either  method  as  before,  we  converge  to  an  extreme  point  on  the  electronic  potential  energy 
hypersurface  which  has  the  aza  nitrogens  defining  the  vertices  of  a  diamond  structure,  see  Figure 
4.  The  pyrrole  nitrogens  are  all  equidistant  from  the  Magnesium  and  the  angles  formed  by  Npy- 
Mg-Npy  are  90°  or  180°,  where  the  Npy  are  the  pyrrole  nitrogens.  This  is  easily  demonstrated 
to  be  an  extreme  point;  displacing  any  one  of  the  pyrrole  nitrogens  0.005A  away  from  the  Mg 
and  reoptimizing  the  structure  leads  to  the  structure  described  in  the  previous  paragraph.  This 
extreme  point  lies  3.0  kcal  per  mole  above  the  true  minima  for  the  INDO/1  structures,  and  some 
3.9  kcal  per  mole  above  the  minima  for  the  ab-initio  structures. 

As  described  above  the  calculated  anion  radical  structures  undergo  a  first  order  Jahn-Teller 
distortion  to  D2h  symmetry.  The  JT  distortion  energy  for  the  INDO/1  structure  was  found  to  be 
2600  cm'1,  or  about  7.4  kcal  per  mole.  This  corresponds  to  the  difference  in  energy  between 
the  “relaxed”  anionic  D2h  structure  and  the  “relaxed”  anionic  D4h  structure.  A  fully  optimized 
anionic  ab-initio  D4h  structure  was  computationally  too  expensive  to  obtain. 

Anion  spectrum.  We  compare  our  computed  spectrum  to  that  obtained  experimentally  by 
Stillman4  et-al  in  Figure  6 b.2S  With  our  choice  of  coordinate  system  the  electronic  ground  state 
is  of  2B2g  symmetry,  from  this  ground  state  the  2Biu  states  are  x-allowed,  the  2AU  states  are 
y-allowed  and  the  2B2u  states  are  z-allowed.  (The  molecular  elongation  is  along  the  y-axis.) 
The  lowest  lying  electronic  excited  state  of  the  anion,  a  2B3g  state,  lies  2900  cm'1  above  the 
ground  state  for  the  INDO/1  structure,  and  4900  cm'1  above  for  the  ab-initio  structure;  it  is 
dipole  forbidden.  In  the  molecular  orbital  picture  of  Figure  3b  this  is  the  OPEN  — *  LUMO 
162g(7r)  [leff(x*)]  — ►  lb^g(ii*)  {leg(ir*)}  excitation,  where  the  bracketed  symmetries  relate  back 
to  the  D4h  parentage  of  these  orbitals.  (This  nomenclature  has  the  HOMO  doubly  occupied, 
and  the  OPEN  orbital  singly  occupied.)  The  effect  of  reducing  the  molecular  symmetry  will, 
of  course,  split  the  Eu  states  in  D4h  into  a  pair  of  non-degenerate  states  in  D2h,  one  being 
x-polarized  and  the  other  y-polarized.  (It  is  obvious  from  Figure  6b  how  the  Q  band  split.)  In 
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the  region  above  790  nm  we  find  two  dipole  allowed  states.  In  the  one-electron  MO  picture 
these  are  characterized  as  promoting  an  electron  occupying  MO  93  into  MO  94,  i.e.  93  — ►  94 
and  the  second  band  is  described  as  93  -  95.  This  pair  of  excitations  would  in  the  neutral 
species  be  the  HOMO  alu(7r)  into  LUMO  eg(-K*)  excitation  that  described  the  Q  band.  In  the 
anjon  these  are  the  HOMO  -+  OPEN  la„(7r)[lalu(7r)]  -»  [le^Tr*)]  and  HOMO  -► 

LUMO  lau(7r)[laiu(7r)]  -*  (7r*)[lep(7r*)J  excitations.  We  agree  with  Stillman4  et-al  that 

t|,c  peaks  observed  at  870  nm  and  790  nm  are  vibrational  overtones  of  the  band  pair  described 
above.  For  both  of  the  optimized  structures  we  compute  the  split  Q  band  of  the  anion  to  be 
(|)  100  low  in  energy  by  approximately  2000  cm'1  and  (2)  the  splitting  of  the  band  to  be  too 
glCal.  Experimentally  the  Q  band  is  interpreted  to  be  split  by  some  300  cm'1  in  the  anion; 

the  INDO/1  structure  gave  a  splitting  of  1400  cm'1  and  the  ab-initio  structure  gave  rise  to  a 
splitting  of  3000  cm'1. 

We  were  curious  as  to  how  the  degree  to  which  the  anion  JT  distortion  affects  the  calculated 
results.  To  investigate  this  we  defined  a  vector  by  the  difference  in  geometry  of  the  D4h  structure 
of  the  neutral  and  the  D2h  structure  of  the  anion  and  proceeded  to  investigate  the  effect  on  the 
calculated  spectrum  as  a  function  of  position  along  this  vector.  We  found  that  the  absolute 
position  of  the  calculated  bands  changed  only  by  a  few  hundred  cm'1,  and  that  the  relative 
positions  of  the  bands  with  respect  to  each  other  also  change  by,  at  most,  a  few  hundred  cm1, 
Wi,h  two  important  exceptions.  The  splitting  of  the  Q  band  was  found  very  sensitive  to  the 
degree  to  which  the  structure  distorted,  and  the  band  we  assign  to  that  observed  at  428  nm  was 
also  quite  sensitive  to  the  degree  of  distortion,  see  Table  IX.  We  found  that  going  from  the 
D4h  structure  to  the  D2h  structure  we  had  to  move  37%  along  the  difference  vector  to  obtain 
a  splitting  of  approximately  300  cm1  for  the  INDO/1  structure,  and  29%  along  the  difference 
vector  to  obtain  the  same  approximate  splitting  for  the  ab-initio  structure.  As  noted,  the  absolute 
posit. otts  of  the  calculated  absorptions  changed  by  a  few  hundred  cm1,  and  for  the  most  part 
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shifted  in  the  direction  that  improved  their  agreement  with  those  experimentally  observed.  Those 
bands  with  parentage  from  the  Eu  bands  of  the  D4h  structure  move  toward  one  another  as  the 
degree  of  distortion  is  lessened,  as  expected.  This  was  most  pronounced  for  the  split  Q  band. 
In  the  one-electron  MO  picture  the  excitations  describing  the  split  Q  band  remain  clean;  that 
is,  there  is  little  mixing  of  other  configurations  into  these  states.  As  we  shall  see  below,  this 
is  not  the  case  for  the  bands  that  lie  higher  in  energy,  see  Table  IX.  The  computed  bands  in 
this  region  occur  at  1360/1274  nm  and  at  1138/1229  nm  for  the  INDO/1  structure,  and  for  the 
ab-initio  structure  they  are  calculated  to  occur  at  1437/1241  nm  and  at  1013/1193  nm.  (Where 
the  first/ second  notation  is  defined  as;  the  first  wavelength  is  for  the  fully  relaxed  D2h  structure 
and  the  second  wavelength  is  for  the  structure  that  best  reproduced  the  Q  band  splitting.)  Overall 
we  find  that  the  split  Q  band  of  the  neutral  shifts  to  lower  energy  in  the  anion  and  is  diminished 
m  intensity,  these  reduction  effects  are  also  observed  in  the  experimentally  obtained  spectrum. 

In  the  region  of  the  observed  spectrum  between  790  nm  to  480  nm  there  is  evidence  of  at 
least  four  bands.  These  occur  with  maxima  at  about  643  nm,  590  nm,  569  nm  and  530  nm, 
from  Figure  6b.  In  this  region,  for  both  optimized  structures,  we  compute  five  dipole  allowed 
bands  with  sufficient  intensity  to  be  easily  observed.  The  relative  intensity  pattern  matches  that 
observed  by  Stillman  et-al,4  Figure  6b.  We  find  two  close  lying  bands  in  the  vicinity  of  the 
shoulder  observed  at  530  nm  in  Figure  6b.  We  feel  it  is  probable  that  the  shoulder  observed 
at  530  nm  in  the  MgPc  spectrum  is  due  to  these  two,  low  intensity,  close  lying  bands.  (The 
observed  spectrum  of  MgPc'  and  ZnPc'  are  remarkably  similar.  We  find  evidence  of  two  bands 
forming  the  shoulder  at  530  nm  in  the  77  K  absorption  spectrum  of  ZnPc‘,41  see  Figure  6c.38 
Given  the  strong  similarities  between  these  two  spectra  we  feel  that  this  lends  credence  to  our 
assignment  of  two  low  intensity  bands  to  this  shoulder.)  The  first  of  the  five  peaks  found  in 
this  region  is  dominated  by  the  94  —  96,  or,  the  OPEN  -►  LUMO+1,  excitation.  This  is  the 
band  we  assign  to  the  peak  observed  at  643  nm  and  we  compute  to  occur  at  812/821  nm  for 


the  INDO/1  structure  and  at  717/763  nm  for  the  ab-initio  structure.  This  is  the  last  absorption 
band  that  we  can  assign  in  this  manner.  All  calculated  bands  lying  higher  in  energy  are  greatly 
mixed  in  the  one-electron  MO  picture.  The  second  band  is  assigned  to  the  shoulder  observed  at 
590  nm  and  computed  at  662/678  nm  (INDO/1)  and  628/662  nm  ( ab-initio ).  The  third  computed 
band  is  attributed  to  the  peak  observed  at  569  nm,  which  we  compute  to  occur  at  532/542  nm 
(INDO/1)  and  at  510/535  nm  ( ab-initio ).  The  fourth  and  possibly  fifth  bands  were  described 
above.  With  this  assignment  the  average  error  is  about  1000  cm'1  for  the  INDO/1  structure  and 
1200  cm 4  for  the  ab-initio  one. 

Experimentally  the  region  below  460  nm  is  complex  with  several  prominences.  Three 
maxima  occur  at  428  nm,  360  nm  and  310  nm.  In  this  energy  range  we  also  compute,  with 
either  structure,  three  comparatively  intense  bands.  For  the  INDO/1  structure  these  occur  at 
407/430  nm,  377/371  nm  and  301/303  nm.  In  the  case  of  the  ab-initio  structure  we  compute 
these  features  to  occur  at  386/426  nm,  358/347  nm  and  281/289  nm.  The  position  of  the  band  we 
assign  to  that  observed  at  428  nm  is  quite  sensitive  to  the  degree  to  which  the  anion  JT  distorts. 
For  both  structures  lessening  the  degree  of  distortion,  i.e.  moving  toward  the  D4h  structure  along 
the  difference  vector,  improves  the  agreement  with  experiment.  From  Table  IX  we  see  that 
when  we  reproduce  the  Q  band  splitting,  the  transition  energy  that  we  associate  with  the  428  nm 
band  is  also  greatly  improved.  Now  if  we  group  the  calculated  bands  by  their  proximity  to  one 
another  and  their  intensity,  albeit  guided  by  the  observed  spectrum,  we  generate  a  pattern  quite 
similar  to  that  observed  experimentally.  In  Tables  VII  and  VIII  these  band  groups  are  designated 
with  {},  and  are  labelled  1-3,  to  indicate  their  correspondence  with  the  observed  spectrum;  we 
also  include  computed  results  down  to  285  nm. 

[PLACE  TABLES  VII,  VIII  AND  IX  HERE.] 

It  is  interesting  to  examine  Tables  VII  -  IX  and  note  how  the  differing  degree  of  distortion, 
principally  of  the  macrocycle,  affects  the  ordering  of  the  states  and  their  intensities.  (A  measure 
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of  how  the  neutral  D4h  species  distorts  upon  reduction,  as  a  function  of  optimization  method,  can 
be  found  in  Tables  I  and  n.)  The  INDO/1  structure  is  “more”  distorted  than  that  of  the  ab-initio 
structure  for  MgPc'.  This  is  easily  seen  from  a  comparison  of  bond  angles  3-10-9  and  7-2-49 
and  bond  lengths  1-6  and  1-19  in  Tables  II  and  I  respectively.  This  has  a  pronounced  effect  on 
the  calculated  spectra  of  the  anion,  especially  for  the  split  Q  band  and  the  computed  428  nm 
peak.  As  is  the  case  for  the  neutral  the  INDO/1  structure  seems  to  produce  a  calculated  spectrum 
for  the  anion  somewhat  in  better  agreement  with  experiment  than  the  ab-initio  structure  does. 

Conclusions 

While  we  could  not  compute  the  Hessian  for  either  molecular  species,  the  computed 
structures  seem  reasonable.  The  neutral  structure  is  assumed  to  be  of  D4h  symmetry;  attempts 
to  warp  the  molecular  plane  or  to  displace  the  Mg  atom  out  of  the  plane  of  the  molecule  lead 
back  to  D4h  structures.  Further,  we  assumed  the  radical  anion  to  be  of  D2h  symmetry,  and  again, 
attempts  to  further  reduce  the  symmetry  of  this  system  lead  back  to  our  initial  D2h  structure. 
Initial  investigation  on  ZnPc  and  ZnPc'  lead  to  structures  in  which  the  Zn  atom  is  displaced 
out  of  the  plane  by  some  0.25  A,  as  is  seen  experimentally.42  At  present  these  systems  were 
studied  at  the  ab-initio  SCF  level  of  theory  using  a  3-21 G  basis  set,  but  our  agreement  with  the 
observed  ZnPc  structure  lend  support  to  the  methodology  used  in  obtaining  the  salient  features 
of  the  MgPc  and  MgPc'  structures. 

Experimentally  there  are  five  strong  bands  observed  between  200  nm  and  800  nm  for  the 
neutral  gas-phase  MgPc  system.35  More  recent  observations  find  that,  at  least,  the  B  band  of  the 
gas-phase  spectrum  is  actually  two  close  lying  intense  bands  split  by  some  3000  cm'1.  We  also 
find  the  B  band  of  the  neutral  species  to  contain  two  close  lying  bands.  The  INDO/1  structure 
reproduces  the  observed  spectrum  somewhat  better  than  the  ab-initio  structure,  perhaps  because 
of  a  larger  Mg-N  distance,  which  is  in  better  agreement  with  experiment.  Additionally  we  find 


a  relatively  intense  band,  not  reported  experimentally,  in  the  L  and  C  region  of  the  spectrum, 
see  Figure  2.  The  lettered  bands  are  described  as  excitations  from  the  ground  state  into  states  of 
E„  electronic  symmetry,  and  are  principally  (tt  ->  n*)  excitations  localized  to  the  macrocycle. 

Our  examination  of  the  neutral  MgPc  spectrum  as  described  above,  does  not  lead  us  to  any 
conclusions  that  differ  from  those  of  References  4,  35  and  36.  We  do  predict  a  weak  z-polarized 
band  hidden  under  the  intense  band  centered  at  350  nm,  see  Tables  V  and  VI,  as  well  as  a  ligand 
to  metal  charge  transfer  band  at  313  nm. 


The  INDO/S-CIS  method  has  been  shown  in  many  examples39, 40, 43-46  to  accurately  repro¬ 
duce  the  low  energy  UV-VIS  spectra  of  compounds  such  as  these.  The  method  consistently 
over  estimates  the  B  band  of  porphyrin  by  some  3000  cm' *  -  5000  cm'^\  this  also  seems  to  be 
true  for  the  phthalocyanines. 


The  calculated  spectrum  of  either  optimized  anion  structure  brings  into  question  the  “good¬ 
ness”  of  the  reported  structure.  The  improvement  of  the  calculated  spectrum,  as  we  move  along 
the  difference  vector  joining  the  D2h  to  the  D4h  structure,  indicates  that  both  optimization  meth¬ 
ods  may  distort  the  molecular  system  to  too  great  a  degree.  The  dramatic  improvement  of  the 
split  Q  band  and  the  band  observed  at  428  nm  as  we  move  from  D2h  back  to  D4h  along  the 
difference  vector  demonstrates  this  claim,  as  well  as  the  general  improvement  of  the  whole  cal¬ 
culated  spectrum.  Even  without  the  “refinement”  of  fitting  the  Q  band  splitting  as  a  function 
of  distortion  to  that  observed  experimentally,  the  calculated  spectrum  as  a  whole  is  good  when 
compared  to  experiment. 


We  find  that  charge  transfer  (CT)  does  not  seem  to  play  an  important  role  below  50000 
cm  1  for  the  anion,  see  Tables  VII  and  VIE.  For  both  species  the  CT  states  we  do  find  are 
macrocycle  to  Mg  charge  transfer  states,  as  would  be  expected  for  Mg+2  Pc'2  and  Mg+2  Pc’3 
molecular  systems.  In  the  case  of  the  anion  there  are  two  states  which  have  a  one-electron 
character  greater  than  40%  Pc  -  Mg(pz),  when  the  experimental  Q  band  splitting  is  reproduced, 
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and  only  one  such  state  is  present  for  the  fully  relaxed  D2h  structure. 

The  four-orbital  model  of  Gouterman  does  describe  the  origins  of  the  Q  and  B  bands  in 
neutral  MgPc.  Unlike  the  porphyrins,  the  Q  and  B  bands  in  Mg  Phthalocyanine  are  not  50/50 
admixtures  of  the  lalu  ->  le*  and  la2«  ->  le*  one-electron  excitations.  The  pseudo  degeneracy 
of  the  loi„  and  la2u  orbitals  found  in  porphyrin  does  not  exist  for  MgPc,  and  this  leads  to  a  more 
intense  Q  band  and  a  less  intense  B  band  for  Pc  as  compared  to  porphyrin.  For  the  anion  the  split 
Q  band  is  described  as  the  laiu  — ►  le*  one-electron  excitation  in  the  neutral;  for  energies  above 
those  of  the  Q  band  the  one-electron  description  of  the  electronic  excitations  loses  clarity,  and 
the  higher  energy  states  are  described  as  admixtures  of  several  such  one-electron  excitations.  For 
this  reason  the  four-orbital  model  is  not  as  complete  a  model  for  the  anion  as  it  is  for  the  neutral. 
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Appendix 


The  Lorentzian  function47  used  in  the  plots  of  Figures  2  and  5  is  given  as 


Vi 


7 rWi 


l  +  4(^f)2 


where  Wi  is  the  half-height  width  of  peak  i,  vs  is  the  sampling  frequency,  i/j  is  the  calculated 
origin  of  peak  i,  fi  is  the  calculated  oscillator  strength  of  peak  i  and  the  summations  i  and  s 
run  over  all  calculated  bands  and  frequency  partitions  respectively. 
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Table  I.  Symmetry  Unique  Bond  Lengths  for  MgPc  and  MgPc  ,  in  Angstroms. 


Atoms3 

MgPc  [D4h] 

6-31G**  INDO/1 

MgPc 

6-3 1G** 

[D2h] 

INDO/1 

1-6 

1.995 

2.098 

1.993 

2.094 

1  -  19 

2.008 

2.109 

3  -  10 

1.317 

1.353 

1.291 

1.339 

3  -  20 

1.356 

1.347 

6  -  10 

1.351 

1.381 

1.362 

1.380 

8  -  9 

1.393 

1.429 

1.389 

1.428 

9  -  10 

1.459 

1.456 

1.476 

1.462 

9  -  14 

1.390 

1.401 

1.386 

1.398 

12  -  13 

1.403 

1.401 

1.397 

1.398 

13  -  14 

1.378 

1.391 

1.383 

1.394 

13  -  17 

1.076 

1.096 

1.077 

1.097 

14  -  18 

1.074 

1.096 

1.075 

1.096 

19-20 

1.351 

1.380 

20  -  21 

1.427 

1.439 

21  -  22 

1.409 

1.439 

21  -  24 

1.406 

1.410 

24  -  25 

1.366 

1.383 

24  -  28 

1.075 

1.097 

25  -  26 

1.420 

1.411 

25  -  29 

1.077 

1.097 

a)  See  Figure  1 . 
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Table  II.  Symmetry  Unique  Bond  Angles  for  MgPc  and  MgPc  ,  in  degrees. 


- 

Atoms3 

«/  ■*  . 

MgPc  [D4h] 

6-3 1G**  INDO/1 

MgPc 

6-31G** 

[D2h] 

INDO/1 

1-6-10 

124.7 

124.2 

124.6 

124.0 

1-19-20 

124.8 

124.1 

3-10-6 

126.6 

126.1 

128.9 

127.1 

3-10-9 

124.1 

126.8 

123.4 

128.7 

3-20-19 

127.1 

126.1 

3-20-21 

124.5 

126.9 

6-1-19 

90.0 

90.0 

90.0 

90.0 

6-1-45 

90.0 

90.0 

6-10-9 

108.3 

107.1 

107.7 

106.6 

7-2-49 

124.6 

128.7 

7-6-10 

110.5 

111.7 

110.9 

112.0 

8-9-10 

106.5 

107.1 

106.8 

107.4 

9-8-11 

121.2 

120.8 

121.2 

120.0 

9-14-13 

117.6 

118.9 

117.8 

119.3 

9-14-18 

120.8 

120.5 

120.6 

120.6 

10-3-20 

125.3 

129.4 

124.6 

128.7 

10-9-14 

132.4 

132.8 

132.0 

132.7 

12  -  13  -  14 

121.2 

121.0 

121.0 

120.8 

14-  13  -  17 

119.7 

119.8 

119.7 

119.8 

19  -  20  -  21 

108.4 

107.0 

20  -  19  -  23 

110.3 

111.8 

20  -  21  -  22 

106.4 

107.1 

20  -  21  -  24 

133.1 

133.4 

21  -  22  -  27 

120.4 

119.5 

21  -  24  -  25 

118.4 

119.7 

21  -  24  -  28 

120.3 

120.0 

24  -  25  -  26 

121.2 

120.8 

24  -  25  -  29 

119.9 

120.2 

a)  See  Figure  1 . 
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Table  nx.  The  Experimentally  Observed,  and  the  Calculated,  Spectrum  of  MgPc. 


Spectrum 

Q 

Ba 

Bi 

Bands  (nm/kK) 

b2 

N 

L 

Vapor*5 

663/15.1 

330/30.3 

— 

— 

280/35.7 

242/41.3 

Solution0 

670/14.9 

— 

361/27.7 

325/30.8 

338/29.6 

281/35.6 

INDO/1 

730/13.7 

— 

325/30.8 

316/32.1 

288/34.7 

239/41.9 

6-3 1G** 

725/13.8 

— 

305/32.8 

302/33.1 

257/34.4 

236/42.0 

a)  The  splitting  of  the  B  band  was  not  observed  in  the  gas-phase  spectrum. 

b)  Reference  35. 

c)  Reference  4. 


C 

219/45.5 

246/40.7 

223/44.9 

221/45.2 
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Table  IV.  Cont. 


44.9  /  1.37 

44  2alu  -+  le* 

11  leg-+lb\u 

n  leg  ^  la*2u 

47.7  /  0.26 

63  leg  ->  1  b*2u 

49.4  /  0.22 

69  la2„  -  2e* 

45.2  /  1.05 

48  2a\ u  — *•  le* 

8  le5  — > 

l^lu 

8  le5  — ► 

l«2u 

48.8  /  0.05 

40  le5  -*• 

16L 

49.9  /  0.08 

26  la2U  - 

18  lo.^ 

12  2eg  — > 

l^u 
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Table  V.  Electronic  Excitation  Spectrum  of  the  Neutral  MgPc  Structure  Optimized  INDO/1 


State 

Symmetry3 

Transition  Energy  (kK) 

Calculated  Observed13 

1 

Ajg 

2,3 

Eu 

13.7  (1.5890)c  Qd 

15.1  [15.9,  16.9]e 

4 

B2g 

23.4 

5 

Aig 

25.6 

6 

A2g 

27.5 

7,8 

Eu 

28.8  (0.0212) 

9,10 

Eg 

29.3 

11 

A2u 

29.9  (0.0331) 

12 

A2g 

30.3  f 

14,15 

Eu 

30.8  (1.1010)  Bi 

30.3  § 

16 

A2g 

31.  lh 

17,18 

Eu 

32.1  (3.3298)  B2 

21 

B2g 

34.7 

24,25 

Eu 

34.7  (0.7004)  N 

35.7 

31,32 

Eu 

36.8  (0.0766) 

38,39 

Eu 

38.1  (0.1286) 

43,44 

Eu 

41.9  (0.5026)  L 

41.3 

47,48 

Eu 

42.9  (0.0432) 

52,53 

Eu 

44.9  (1.3680)  C 

45.5 

a)  (x,y)  transforms  as  eu  and  (z)  transforms  as  a2u  in  EUh- 

b)  From  Reference  35. 

c)  Values  in  parenthesis  are  calculated  Oscillator  Strengths. 

d)  Spectroscopic  labels  for  the  degenerate  transitions,  Q,  Bi,  B2,  N,  L  and  C. 

See  Reference  35. 

e)  Vibrational  fine  structure,  see  text,  as  well  as  References  4  and  35. 

f)  Byond  30300  cm*1  only  states  with  calculated  oscillator  strength  are  reported. 

g)  The  absorption  is  broad  ,  see  Reference  4,  and  is  assigned  to  calculated  states  (7-18); 
see  also  Figure  2a. 

h)  The  only  state  below  50000  cm'1  with  significant  charge  transfer  character; 

50%  Pc  — ►  Mg  (7 r  — ►  7 r*). 
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Table  VI.  Electronic  Excitation  Spectrum  of  the  Neutral  MgPc  Structure  Optimized  6-3 1G** 


Transition  Energy  (kK) 

State 

Symmetry3 

Calculated 

Observed1* 

1 

Aig 

2,3 

Eu 

13.8  (1.5504)c  Qd 

15.1  [15.9,  16.9]e 

4 

B2g 

23.7 

5 

Aig 

25.6 

6 

A2g 

27.3 

7,8 

Eu 

28.6  (0.0374) 

9 

A2g 

30.9  f 

13 

A2u 

31.9  8  (0.0364) 

15,16 

Eu 

32.8  (1.2856)  Bi 

30. 3h 

18,19 

Eu 

33.1  (1.6932)  B2 

22,23 

Eu 

34.4  (1.8004)  N 

35.7 

29,30 

Eu 

36.7  (0.0038) 

38,39 

Eu 

38.9  (0.2190) 

43,44 

Eu 

42.0  (0.1999)  L 

41.3 

45,46 

Eu 

42.8  (0.0594) 

51,52 

Eu 

45.2  (1.0470)  C 

45.5 

a)  (x,y)  transforms  as  eu  and  (z)  transforms  as  <2211  in  D^. 

b)  From  Reference  35. 

c)  Values  in  parenthesis  are  calculated  Oscillator  Strengths. 

d)  Spectroscopic  labels  for  the  degenerate  transitions,  Q,  Bi,  B2,  N,  L  and  C. 

See  Reference  35. 

e)  Vibrational  fine  structure,  see  text,  as  well  as  References  4  and  35. 

f)  The  only  state  below  50000  cm'1  with  significant  charge  transfer  character; 
50%  Pc  — ►  Mg  (7T  — *■  7T*). 

g)  Byond  31900  cm'1  only  states  with  calculated  oscillator  strength  are  reported. 

h)  The  observed  band  is  broad,  see  Reference  4,  and  is  assigned  to  states  (7-19); 
see  also  Figure  2b. 


Table  VII.  Electronic  Excitation  Spectrum  of  the  Anion  MgPc  Structure  Optimized  INDO/1 


State 

Symmetry3 

Transition  Energy  (kK) 

Calculated  Observedb 

1 

®2g 

2 

®3g 

2.89 

3 

Au 

7.36  (0.0821)c 

10.5 

4 

Biu 

8.79  (0.0063) 

10.8  [11.5,  12.7]d 

5 

Au 

12.3  (0.4637) 

15.6 

6 

B3g 

14.6 

7 

Biu 

15.1  (0.2313) 

16.9 

8 

B2g 

18.5 

9 

Biu 

18.8  (0.4570) 

17.6 

10 

B3g 

19.9 

11 

Au 

20.3  (0.1116) 

18.9 

12 

Biu 

21.2  e  (0.2110) 

14 

Au 

22.7  (0.0066) 

17 

Biu 

23.5  (0.0143) 

19 

Biu 

24.5  (0.1454) 

23.4  {l}f 

21 

Au 

26.2  (0.0464) 

22 

Biu 

26.5  (0.4601) 

27.8  {2} 

25 

b3u 

27.2  (0.0003) 

29 

Biu 

28.4  «  (0.0999) 

30 

Au 

28.4  (0.0522) 

31 

Biu 

28.6  (0.0018) 

34 

Au 

29.5  (0.1246) 

35 

Biu 

29.7  (0.1171) 

37 

Au 

30.5  (0.2676) 

41 

b3u 

31.8  (0.0316) 

44 

Biu 

33.1  (0.0008) 

46 

Au 

33.2  (0.8541) 

32.3  {3} 

a)  With  our  choice  of  coordinates  Biu  is  (x)  allowed,  Au  is  (y)  allowed  and  B3u  is  (z)  allowed. 

b)  From  Reference  4. 

c)  Values  in  parenthesis  are  calculated  oscillator  strengths. 

d)  See  text.  These  are  assigned  as  overtones  of  the  transition  observed  at  10500  cm'1. 

e)  Beyond  21200  cm'1  only  states  with  computed  oscillator  strength  are  reported. 

f)  See  text.  Each  observed  band  is  assigned  to  several  calculated.  See  also  Figure  5a. 

g)  The  only  state  below  50000  cm'1  with  significant  charge  transfer  character; 

48%  PC  ->  Mg  (7T  7T*). 
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Table  VIIL  Electronic  Excitation  Spectrum  of  the  Anion  MgPc  Structure  Optimized  6-3 1G** 


Transition  Energy  (kK) 


State 

Symmetry3 

Calculated 

Observed6 

1 

®2g 

2 

®3g 

4.93 

3 

Au 

6.96  (0.1231)' 

10.5 

4 

Biu 

9.87  (0.0037) 

10.8  [11.5,  12.7]d 

5 

Au 

13.9  (0.4463) 

15.6 

6 

B3g 

15.4 

7 

Biu 

15.9  (0.2255) 

16.9 

8 

B2g 

19.4 

9 

Biu 

19.6  (0.3652) 

17.6 

10 

Bjg 

20.6 

11 

Au 

20.9  (0.1043) 

18.9 

12 

Biu 

21.9  c  (0.2072) 

14 

Au 

23.2  (0.0014) 

17 

Biu 

24.0  (0.0381) 

20 

Biu 

25.9  (0.1312) 

23.4  [1  }f 

21 

Au 

26.3  (0.0319) 

22 

B3g 

27.1 

23 

Biu 

27.6  (0.2831) 

27.8  {2} 

25 

Biu 

28.2  (0.1162) 

27 

Au 

29.1  (0.0025) 

29 

Biu 

29.9  g  (0.0678) 

30 

Au 

30.1  (0.0548) 

33 

Biu 

30.9  (0.1409) 

36 

Au 

32.0  (0.2151) 

32.3 

40 

Biu 

33.5  (0.0562) 

41 

b3u 

34.0  (0.0026) 

42 

Biu 

34.1  (0.0482) 

43 

Au 

34.7  (0.0137) 

47 

B3u 

35.2  (0.0001) 

48 

Au 

35.4  (0.8878) 

{3} 

49 

Au 

35.8  (0.2839) 

a)  With  our  choice  of  coordinates  Bi„  is  (x)  allowed,  Au  is  (y)  allowed  and  B3u  is  (z)  allowed 

b)  From  Reference  4. 

c)  Values  in  parenthesis  are  calculated  oscillator  strengths. 

d)  See  text.  These  are  assigned  as  overtones  of  the  transition  observed  at  10500  cm'1. 

e)  Beyond  21900  cm  1  only  states  with  calculates  oscillator  strength  are  reported. 

f)  See  text.  Each  observed  band  is  assigned  to  several  calculated.  See  also  Figure  5b. 

g)  The  only  state  below  50000  cm*1  with  significant  charge  transfer  character;  48%  Pc  -+Mg  (tt  -►  tt*). 


26 


Table  IX.  The  Allowed  Prominent  Bands  of  MgPc'  as  a  Function  of  the  %JT  Distortion  and 


Optimization  Method.  (Energy  in  units  of  nm  /  kK.) 


INDO/1  structure 

ab-initio  structure 

D2h 

37%  D2h 

D2h 

29%  D2h 

Experiment 

1359  /  7.36 

1274  /  7.85 

1437  /  6.96 

1241  /  8.06 

952  /  10.5 

1138  /  8.79 

1229  /  8.14 

1013/9.87 

1193  /  8.38 

923  /  10.8 

813  /  12.3 

820  /  12.2 

719/  13.9 

763  /  13.1  ' 

643  /  15.6 

662  /  15.1 

680  /  14.7 

629  /  15.9 

662/  15.1 

592  /  16.9 

529  /  18.9 

543  /  18.4 

510/  19.6 

535  /  18.7 

569  /  17.6 

493  /  20.3 

490  /  20.4 

478  /  20.9 

483  /  20.7 

529  /  18.9 

472/21.2 

467  /  21.4 

457  /  21.9 

459/  21.8 

529  /  18.9 

408  /  24.5 

429  /  23.3 

386  /  25.9 

426  /  23.5 

428  /  23.4 

377  /  26.5 

370  /  27.0 

362  /  27.6 

347  /  28.8 

360  /  27.8 

301  /  33.2 

303  /  33.0 

282  /  35.4 

289  /  34.6 

310/32.3 
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Figure  Captions 


Figure  1.  Magnesium  Phthalocyanine,  and  the  atomic  numbering  scheme  used  for  the 
Tables  and  text.  The  elongation  for  the  anion  is  along  the  y-axis,  which  is  defined  by  atoms  6 
and  32.  The  z-axis  is  perpendicular  to  the  molecular  plane. 

Figure  2.  The  calculated  stick  spectrum  of  neutral  MgPc.  a)  For  the  INDO/1  optimized 
structure  and  b)  for  the  6-3 1G**  optimized  structure.  The  lines  are  broadened  using  a  Lorentzian 
line  shape  and  a  width  at  half  height  of  250  cm'1,  see  Appendix. 

Figure  3.  The  MO  diagrams  for  MgPc  and  MgPc',  a)  and  b)  respectively.  The  right  most 
column  is  the  orbital  energy  in  eV,  the  ( M  N)  label  the  MO  and  its  symmetry  in  the  D2h  point 
group,  respectively.  The  symmetry  labels  N  run  from  1-8  and  are  respectively  Ag,  Au,  Big, 
Biu,  B2g,  B2u,  B3U  and  B3g,  i.e.  (12  8)  would  be  MO  M=12  of  symmetry  B3g.  The  right  most 
columns  label  the  orbitals  as  of  either  sigma  or  pi  type. 

Figure  4.  A  symbolic  representation  of  the  possible  Jahn-Teller  distortions  of  MgPc'.  Na 
represents  the  aza-nitrogens  and  Np  the  pyrrole  nitrogens.  The  square  of  Np,  the  top  and  bottom 
figures,  are  calculated  to  be  the  most  stable.  The  four  Np  then  form  two  possible  rectangles. 
The  other  two  structures,  the  left  and  right  most  figures,  with  the  four  Na  diamond  like,  and  the 
four  Np  square  like,  represent  saddle  points  connecting  the  top  and  bottom  structures. 

Figure  5.  The  calculated  spectrum  of  MgPc  anion,  a)  For  the  ENDO/1  optimized  structure 
and  b )  for  the  6-3 1G**  structure.  The  lines  are  broadened  using  a  Lorentzian  line  shape  and  a 
width  at  half  height  of  250  cm'1,  see  Appendix. 


Figure  6.  The  observed  spectra  of  MgPc  a),  MgPc'  b)  and  ZnPc'  c). 
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